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Abstract 

In this paper we complete a description of calculation of the one-loop amplitude for 
e + e~ — * ff process started in CERN-TH/2001-308. This study is performed within the 
framework of the project CalcPHEP. Here we add QED subsets of the one-loop diagrams 
and the soft-photon contribution. The formulae we derived are realized in two indepen- 
dent FORTRAN codes, eef f Lib, which was written in an old fashioned way, i.e. manually, 
and another one, created automatically with an aid of s2n_f (symbols to numbers) soft- 
ware — a part of CalcPHEP system. We present a comprehensive comparison between 
the two our codes as well as with the results existing in the world literature. 
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Introduction 



Recently, detail reports on process e + e~ — > tt — > 6/ become an active subject for energy of 
future electron linear colliders. This process will be one of main process and therefore must 
be theoretically studied profoundly (see for example the review [1]). 

In this connection we consider a new calculation of e + e~ — > // process at the one-loop 
level made with an aid of computer system CalcPHEP, where all the calculations from the 
Lagrangians up to numbers are going to be eventually automatized, (see [2]). 

Electroweak (EW) parts have been calculated in [3] and a very good agreement with 
FeynArts [4] and [5] were found. 

In this paper we added lacking in [3] QED corrections. Our strategy in the descriptions 
of the QED part is the same as in our first paper; many definitions and notations from it are 
used here. References to an equation of the first part will be denoted as (I.S.eq) with S and 
eq being Section and equation numbers of Ref. [3], correspondingly. 

This paper is organized in a similar fashion as [3]. 

In Section 1, we briefly remind the structure of one-loop amplitudes 

Section 2 contains explicit expressions for all the QED building blocks which were not 
covered in [3]: QED vertices, AA and ZA boxes. 

Section 3 contains the total scalar form factors of the one-loop amplitudes, now with all 
QED additions. 

In Section 4 we present explicit expressions for helicity amplitudes made of total scalar 
form factors at one-loop level. 

Section 5 is an Annex containing some additional expression for different QED contri- 
butions that might be derived analytically. They are not in the main stream of our paper: 
Lagrangian — > scalar form factors — > helicity amplitudes — * one-loop differential cross-section. 
However, they are useful for pedagogical reasons, and their coding in complimentary FORTRAN 
branches of eef f Lib provided us with powerful internal cross-checks of our codes for numerical 
calculations. Actually, eef f Lib version of February'2002 has three QED branches. 

Finally, Section 6 is a revised version of Section 5 of [3] in which we present again results 
of a comprehensive numerical comparison between eef f Lib and ZFITTER. The reason for this 
revision is due to debugging of the December'2001 version of eef f Lib resulting in a little 
change of our numbers beginning 4th or 5th digits. In this paper we also present a comparison 
with our another code, which was created automatically using s2n_f software. We also present 
a comprehensive comparison between the results derived with two our codes and the results 
existing in the world literature. In particular, we found a high precision agreement with 
FeynArts results up to 11 digits for the differential cross-sections with virtual corrections, and 
with resent results of [5] within 7-8 digits even with soft photons included, see [6]. 
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1 Amplitudes 



We work in the LQD basis, and the final-state fermion masses are not ignored as in previous [3] . 
The electron mass is ignored everywhere, but arguments of logs. Also we work in the 
gauge. We checked the cancellation of ^-dependent terms in three gauge-invariant subsets of 
diagrams separately. The first subset is the so-called cluster in the QED sector (or A cluster, 
see definitions below), the second and third are AA boxes and ZA boxes, correspondingly. 
In the LQD basis, the 7 and Z exchange one-loop amplitudes have the following structure: 

AnQ e Qf 



A™ A = i^^a(s) llx ®^, (1.1) 



and 



+/ e (3) ^7+ ® l„F LQ (s, t) + 5 e 5 t l„ ® 1„F QQ (s, t) 

+I^I^ lfll+ ® i-imtD^) F LD (s,t) + 8 e I?\ ® (-im t D^) F QD (s,t)|, (1.2) 

where untilded and tilded form factors are related by Eqs. (1. 1.11). Like Part I, we present 
all the explicit expressions in term of untilded quantities. Furthermore, 

a(s) = Ty— (1.3) 

i-^l*)-*] 

is the fermionic component of the running QED coupling a(s) and 

Xz{s) = r~ r~r~ (L4) 

M z 

is the Z/7 propagator ratio with an s-dependent (or constant) Z width. 



2 Building Blocks (QED part) 

2.1 The Zff and 7// vertices 

First of all we have to add vertex QED building blocks to the scalar form factors of Eq. (1.2.60) 
and finally to the complete scalar form factors of Eqs. (1. 3. 117). 

The total vertex scalar form factors 7W and Ztt Eqs. (1.2.60) are now sums over all 
bosonic contributions B = A, Z, W, H, since we add the diagram with virtual 7 = A. 

All the 24 components of the total form factors in the LQD basis look like: 

(*) = (*) + (*) + *ZST (*) + F2$j (s), (2.1) 

The A cluster was formed using the same philosophy as in [3], see Eqs. (I.2.54)-(I.2.59). 

Note, that F^ A (s) and F2 H (s) are equal to zero. 
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2.1.1 Library of QED Form Factors for Att clusters 

Up to one- loop level, there are two diagrams, which contribute to the A cluster, see Fig. 1. 




+ 





Figure 1: A cluster. One fermionic self-energy diagram in brackets gives rise to the counter 
term contribution depicted by the solid cross. 

Since after wave function renormalization, the scalar form factors became UV-finite, instead 
of Eq. (1.2.61), we have for all 6 form factors which are also separately gauge-invariant: 



F~i(z)A _ <py(z)A 
1 I J I 1 

where I = L,Q, D. Individual components are: 
TT = 0, 

?Z A = Qt s w{ 2 ( s - 2m t) Cb( - m t , ~ m t , -s; m u o, m) 

-3Bq( - s;m t ,m t ) +SBq( -m 2 ;m t ,0) - 4m 2 B 0p ( -m 2 ;0. ///,) \. 
Bq( - s; m t , m t ) - £jf ( - m 2 ; m t , 0) 

Bq\ - s; m t , m t ) - B%( - m 2 ; m t , 0) 

Bq( ~ s; m t , m t ) - Bfi( - m 2 ; m t , 0) 



D 


Qt s w ^ 




- s; 




spzA 

L 


= J~Q + Qt s w 


8m 2 
A 3r 


<( 


<rzA 

Q 


T ^A _ n 1 2 


8m 2 


r (3) r 


A 3r 


^ L 


*pzA 
D 


Q 2 t*l 2v t 

4 3) ±* 




- s; 



( - s; mt, mt) - B (-m t ; m t , 0) 



with 



A 3r = 4m? 



s . 



2.1.2 Scalar form factor for electron case 

Aee cluster is described by only one scalar form factor: 



F A ' e (s) 



2„2 



Qe S 



2sC ( - m 2 e , -m\, -s; m e , 0, m e ) 



-3Bq( - s; m e , m e ) + 3Bq( - m\; m e , 0) - 4m 2 e B 0p ( - m 2 e ; 0, m e ) 



(2.2) 



(2-3) 



(2.4) 



(2.5) 
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2.2 Amplitudes of QED boxes 



The contributions of QED AA and ZA boxes form gauge-invariant and UV finite subsets. In 
terms of six structures (L, R) ® (L, R, D) they read: 



\ d+c 
gAA(ZA) \ _ 



iAA(ZA)9_ 

"'nnvm 



[ 7m7+ ® ?tt' ZA) (s, t, u) + [ 7 ^7 + ® 7,7-] ^ A(ZA) (a, * , «) 



+ [1,1- ® 7^7+] ^" (ZA) («, t, u) + [ 7 „ 7 _ ® 1,1 A Hi {ZA) {s, t, u) 



(2.6) 



+ [ 7m 7+ ® (-im*/^)] TlT A) (s, t, u) + [ 7/l7 - ® Hm t ID,)} T A R T A) («, t, u) 



where for shortening of presentation we factorize out normalization factors: 



UAA _ 4 /O2.o2 
^norm ~~ WeVt j 



^,za _ s wQeQt 



(2.7) 



For completeness and subsequent use we remind k^ rm appearing in Eq. (1.2.95): 

1 



l.zz __ 

norm 32cl 



{2.1 



2.2.1 AA-box contribution 

There are only two AA diagrams, direct and crossed: 



■:,+ 



A 



t 



WWVWWWWWf 



+ 




— *- — AAAAAAAAAAAA/Wf — ► — 

e A t e A t 

Figure 2: Direct and crossed AA boxes. 

The six form factors of AA boxes might be expressed in terms of only four auxiliary 
functions T\ and H,\ 2 3: 



F^( s ,t,u) = T^(s,t,u) = n 1 (s,t)-n 1 (s,u) + n 2 (s,t) + n 3 (s,u), 

F%(s,t,u) = F^(s,t,u) = Hi(s,t)-H 1 (s,u)-H 2 (s,u)-H 3 (s,t), 
F%(s,t,u) = F%(8,t,u)=F 1 (s,t)-r 1 {s,u). (2.9) 



The auxiliary functions are rather short: 
1 s [ 1 / 



2 I ^4r 



t_J AA (-s, -t;m e ,m t ) 



+ts 



+2- 



*3r 



s C ( - ml, —ml, -s; 0, m e , 0) + (s — 2mfj C ( - m 2 , -m 2 , -s; 0, m t , 0) 
2m 2 C ( - m 2 t , -m 2 t , -s; 0, m t , 0) + P F ( - m 2 t ; m u 0) - B%{ - s; 0, 0) 



t 

+2 r 



Bq{ - t; m e , m t ) - Bq( - m\; m t , 0) 



(2.10) 



Hi(s,t) = -t_ -J AA (-s,-t;m e ,m t ) - C ( 



m e , — m 



! , -t;m e ,0,m t ) 



4A 



4r 



Mr 



^aa(-s, -t;m e ,m t ) 



+2m t 2 ll-2- 



Bq( ~ t; m e , m t ) - Bq( - m 2 ; m t , 0) 



(2.11) 



n 2 (s,t) 



4A 4r 

xC ( - m\, -m 2 , -s; 0, m t , 0) 



- 2m 2 s + (s - Am*) (s + 2t_ - (st+ + 2tt_) 



Mr 



+ (s + 2i_) - *+^~) s Co( - "ie» ~ m e ? -s; 0, m e , 0) 



-2*. 



-4m 2 



Here 



4A 



■4r 



P F (-t;m e ,m,)-P F (-s;0,0) 
Bjf ( - t; m e , m t ) - - m 2 ; m t , 0) 
(s + 2tJ) t 3 _J AA (-s, -t; m e , m t ) . 

A 4r = — tu + mf , 



(2.12) 
(2.13) 

(2.14) 



and J AA {Q 2 ,P 2 ;M 1 ,M 2 ) is due to a procedure of disentengling of the infrared divergences 
from D . Its explicit expression reads (P 2 > 0, Q 2 < 0, and M 1 is ignored everywhere but In): 



(P 2 + M 2 f ( P 2 \ 1 



In 2 



M 2 



+ >4 + M)- 2 n 2 



P 5 



P 2 + M 2 



+ in In 



-<2 2 , 

(P 2 + M 



-In 5 



,2-\ 



Ml Ml 



-Q 2 ) 

(2.15) 



Moreover, the relevant infrared divergent C function (P 2 > again), is 



C^(-M 2 1 ,-M 2 ,P 2 ;M 1 ,\,M 2 ) = ^Jj 



+ M 2 ) 



In 



(P 2 + M 2 
M\Ml 



P 2 

In — - 2Li 2 



P 2 + Mo 2 
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2.2.2 ZA box contribution 



In i?£ gauge there are eight ZA boxes, however, since electron mass is ignored, only four 
diagrams without O contribute: 

Z i „+ A * 



o + 



wwvwvw* 



-*2 — AAAAAAAAAAAA/Wf — 

e A t 



■:,+ 







t 



— AAAAAAAAAAAA/Wf — +- 

e A t 

7 t 




+ 



+ 



+ 



wwwwwwww 



— •V/WWWWv — 

e Z t 



t 



wwwwwwww 



0° 

A 




z 

A 



+ 



At e 
Figure 3: Direct and crossed boxes. 




The six relevant scalar form factors are conveniently presentable in form of differences of 
t and u dependent functions: 



(2.17) 
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where index I J is any pair of L,R © L, R, D. The 12 J~ff functions depend on 6 auxiliary 
functions by means of equations where the coupling constants are factored out: 



-TZZA 

LL 


(s,t) 


= Ve^tGl + <V*t& (S,t) , 


^ (*, 


■u) 




, U) + (T £ (T t H 2 (s, U) . 


J~R.Il 


(s,t) 


= S e 5 t Q! (s, t) + 5 e a t Q 2 (s, t) , 


RR ( S - 


■u) 


= S e a t Hi(s 


, u) + 5 e 5 t H 2 (s 


u) , 


qzzA 

LR 


(s,t) 




^ (*, 


■u) 


= (r e S t Gi(s, 


u) + o e a t Q 2 (s 


u) , 


rzA 
RL 


(s,t) 


= <5 e <W (s, t) + 5 e a t H 2 (s, t) , 




u) 


= ^> e o t Qx (s, 


u) + 8 e 8 t g 2 (s, 


«) , 


rzA 
LD 


(s,t) 


= Ve^iFl ( S > + <7e S tF- ( S ' > 


J~* LD ( S - 


u) 


= 






qzZA 
RD 


(s,t) 


= ^A-? 7 ! (s, i) + 5 e ^^2 (s, , 




u) 


= S^J 7 ! (S, 


u) + 5 e 5 t< F 2 (s, 


u). 



Finally, we present these 6 auxiliary functions: 



F 2 (s,t) 
Hi (s,t) 



1 s 



8 A 



4r 



f- 



-Rz + 2 ) J ZA (-s, -t; m e , m t ) 

s 



-4C ( - to?, TOg, -t; TO t , M z , to 6 ) 
s_t_ r 

+2— — tJ ZA (-s, -t; to 6 , TOt) + 2tC ( - to?, to?, to*, M z , m e ) 
A 4r L 

-t_C ( - to?, -to?, -s; M z , TO e , 0) - t+C ( - ™?, -to?, -s; M z , to 4 , 0) 
-s + C ( - TO?, -to?, -s; M z , TO e , 0) - s_C ( - to?, -to?, -s; M z , TO t , 0) 



-2t C ( - TO?, -TO?, -t; TO t , M z , TO e ) - 2-^ 

6 — 



-25 F ( - t; TO e , TO t ) + 5 F ( - to?; M z , to 4 ) + 5 F ( - to?; m t , 0) 



M Z C ( - m t , TO e , -t; mt, M z , to 6 ) 



*3r 



(M? - 4to?) C ( - to?, -to?, -s; M z , to 4 , 0) 



+25 F ( - s; M z , 0) - <( - to?; M z , m t ) - <( - to?; m t , 0) 



1 s 



8A 



4r 



STO7 



Jza{s, ~t; TO e , TOt) - 2tC ( - TO?, -TO?, TOt, M z , TO e ) 

M 2 

Jza(s, -t, TO e , TOt) - 2— ^C ( - TO?, -TO?, -t; TOt, M z , TO e ) 

1/ — 



(2.18) 
(2.19) 



+ 



Mr 



Jza(~S, -t] TO e , TOt) + 2C ( - TO?, -TO?, -t; TOt, M z , TO e ) 



+4to?C ( - to?, -to?, -s; M z , TOt, 0) + t- \ R Z - 2 + 

,2 _2 _ S/ r _ \ o^/ ™2 _2 



A4 r 



^za(-s, -t;m e ,m t ) 



+2C ( - to 4 , TO e , -t; TO t , M z , TO e ) - 2C ( - TO t , -m t , -s; M z , to 4 , 0) 
t 



+s_ (s + 2t_) 



Mr 



C ( - TO t , -TO t , -s; M z ,TO t ,0) 
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-C ( - m 2 , -m 2 e , -s; M z , m e , 0) 

+ (R z - 1) (s + 2f_) C ( - m ( 2 , -m 2 t , -s; M z , m t , 0) 
-2Bg( - s; M z , 0) + 2B£( - t; m e , m t ) 



s + C ( - m 2 e , -m 2 e , -s; M z , m e , 0) 



(2.20) 



+2— 

^ 



W 2 (M) = s| 



2<( - t; m e , m«) - <( - m 2 t ; M z , m t ) - B*( - m 2 ; m t , 0) 

Jza(s, -t; m e , m t ) - C ( - m 2 e , -m\, -t; m e , 0, m t ) 

mf (\ 



4s_ 

+C ( - m 2 , -m 2 , -t; m t , M z , m e ) 



mi ( 1 T . . 
+ — "i e , "^J 

8 \s 



m 



2 r 



Mr 



£?i (s,*) = s 



t- 



Jza(s, -t; m e , m t ) + 2C ( - m 2 , -m 2 e , -t; m t , M z , m e ) 
- Jza(s, -t; m e , m t ) + C ( - m 2 , -m 2 , -i; m e , 0, m t ) 



(2.21) 



4s. 

-C ( - m 2 , -m 2 , -t; m t , M z , m e ) 



m; 



- ^JzAis, -f ; m e , mt) 
8s 



8A 4r 



Jza{s, -t; m e , m t ) + 2C ( - m\, -m 2 e , -t; m t , M z , m e ) 
Jza(s, -t;m e ,m t ) 



+ mt 



mfM 2 
2t, — z - 



-At + C ( - ml, -m 2 e , -t; m t , M z , m e ) + 2m 2 C ( - m 2 t , -m 2 t , -s; M z , m t , 0) 



S-t-t A 



\ 4r 



tJ ZA (-s, -t; m e , m t ) + 2tC { - mf, -m 2 e , -t; m t , M z , m e ) 



-t+C ( - ml, -ml, -s; M z , m t , 0) - t_C ( - m 2 e , -m 2 e , -s; M z , m e , 0) 

+t s+C ( - m 2 e , -m 2 e , -s; M z , m e , 0) + (s+ - 4m 2 ) C ( - m 2 , -m 2 , -s; M z , m t , 0) 

+2B%( - s; M z , 0) - 2B$( - t; m e , m t ) 
m 2 t 



+2- 



t. 



M Z C ( - m t , -m e , -t; m t , M z , m e ) 



-2B*( - f, m e , m t ) + B*{ - m 2 ; M z , m t ) + B£( - m 2 ; m t , 0) 



smt 



t 



A 



Ar 



-Jza{~s, -t;m e ,m t ) 
is 

J ZA {-s, -t; m e , m t ) + 2C ( - m 2 t , -m 2 , -t; m t , M z , m e ) 



where new notation were introduced for invariants 



s± = s±M 2 , t± = t±ml , 



(2.22) 

(2.23) 
(2.24) 
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and for the new functions J IJ (Q 2 , P 2 ; M u M 2 ) (an analog of J AA (Q 2 , P 2 ; M u M 2 ) Eq. (2.15)): 



J AZ (Q 2 ,P 2 ;M 1 ,M 2 ) 



1 



P 2 + Mf 



M 2 M 2 



(P 2 + M 2 ) 



2^2 



v 2 -25) 



2.2.3 Box— Born interferences 



Any box, describing by the amplitude Eq. (2.6), interfering with 7 and Z exchange tree level 
amplitudes, gives rise to two contributions to the differential cross-sections, which are useful 
for internal cross-checks: 



Cr BOX(g)BORN 7 0^ 8QeQ*R' e '! 



(s + t_) 2 + sm\ 



{T LL + T R 



+ (sm 2 t + t 2 _) {T LR + T RL ) - 2m 2 (st + t 2 _) {T LD + T RD ) ) J, (2.26) 



Cboxixibornz SRe^ ^ 



(s + tS) 2 + sm\ 



5 t {a e F LL + 5 e F RR ) (2.27) 
+2 (s + t_) 2 a t a e T LL + 2t 2 _a t 6 e F RL + (sm 2 t + *i) <S t [a e T LR + 5 e ^ M ) 
+2sm 2 a t {o e T LR + 5/^) - 2m 2 (st + t 2 ^ v t {o e T LD + b e T RD ) )xz*j- 

3 Total scalar form factors of the one-loop amplitude 

Adding all contributions together, we observe the cancellation of all poles. The ultraviolet- 
finite results for six scalar form factors, replacing EW result Eq. (1.3.118), are: 



F LL (s,t,u) 
Fql (s,t,u) 
F LQ (s,t,u) 
F QQ ( s ,t, u ) 

F LD (s, t, u) 

F QD (s, t, u) 
where 



TT (s) + T A,e (s) ] + r tt (s) + ^ (s) + 16fcF «,x (S; tj u) j 

^ ee (*) + ^ + ( s ) + k H tt (s) + T* L (s) + 16k T B Q L X (s, t, u) , 



Tf e (s) + T 



■A,e 



+ J*< (s) + k Tr (s) + J% (s) + 16kT°° x (s, t, u) , 



Tf (s) + k FT (s) + 16kF°° x (s, t, u) , 



(3.1) 
(3.2) 



k = c 2 w (R z - 1) . 
For I J = LL component of box contribution one has: 

Tf° x (s, t, u) = k AA Tf A (s, t, u) + k ZA T ZA (s, t, u) + k zz T z f (s, t, u) + k ww T™ w (s, t, u) (3.3) 
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and for the other components IJ = LQ, QL, QQ, LD, QD of box form factors the WW box 
does not contribute. Moreover, 

r£$S(*)= £ *7S£(«), (3-4) 

B=A,Z,H,W 

except J=l A (s) = and F} H (s) = 0. 



4 Process eett in the helicity amplitudes 

According to the analysis of the EW part in [3] and presentation of the QED part here, we 
have the complete answer for the amplitude of our process. 

The aim of this section is to adapt the helicity amplitude technigues for the description 
of our process. We produced an alternative analityc answer for the same amplitude using the 
method suggested by Vega and Wudka (VW) [7]. 

In general, there are 16 helicity amplitude for any 2f — > 2f process. For the unpolar- 
ized case and when the electron mass is ignored, we are left with six independent helicity 
amplitudes, which depend on kinematical variables and our six form factors: 



A++++ 
A+-+- 

A+—+ 

A+— 

A- +++ 



A- 



A. 



•++- 



-+ 



Here 



0, A +++ - = 0, 

s (1 - COS??) (Q e Q t F G G + XzS, 
S (1 + COS??) (Q e Q t F GG + XzS, 



A++-+ = 0, A ++ — = 0, 
(l + P t )4 3) F QL + S t F ( 

[l-P t )li 3) F QL + S t F QQ 



QQ 



2 y r sm t sin ■& [Q e QtF G G + Xz$e 
A-+— = -2^m t sin $(Q e Q t F GG 



4 3) F QL + 5 t F QQ + ±s[3 t 2 li 3) F QD 



+X> 



2I^4 3) F LL + 2I^S t F 



t 1 LQ 



S e 4 3) F OL + S e S t F 



e^t 1 QQ 



+^A 2 4 3) (2I^F LD + 5 e F QD ) 
= s(l + costf) (Q e QtF GG 



+x 



(1 + fa) (2lV>I®F LL + 6 e ll 3) F QL ) + 5 t (2I^F LQ + 5 e F ( 

S (1 - COS$) (QeQtF GG 



+Xz 

0, 



(1 - p t ) 4 3) (2I^F LL + 5 e F QL ) + S t (2iPF LQ + 5 e F 
A—+- = 0, A—+ = 0, A = 0. 



<>Q 



cos$ 



t-mf + 



27 sp t 



(4.1) 



(4.2) 
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and for the amplitude A\ i \ j \ k \ l each index X^j^i) takes two values ± meaning twice projection 
of spins e + ,e~,t,i onto their corresponding momentum. The differential cross-section for the 
unpolarized case is: 

PtN c ]T ^W. • (4.3) 



dcosfi s 3 . , 

We checked, that this expression is analytically identical to Eq. (1.4.122). The expression 
Eq. (4.3) contains, however, spurious contributions of the two- loop order (squares of one-loop 
terms), which one should supress, since we would like to have a complete one-loop result. 

This may be achieved with a simple trick. First of all let us note, that if all form factors 
are: F I3 = 1 for I J = LL, LQ, QL, QQ and F I3 = for I J = LD, QD, we have the tree 
level. At the one-loop level LL, LQ, QL, QQ form factors may be represented as: 



and 



F " = ik F "' (45) 



for I J = LD, QD. 

Instead of Eq. (4.4) for the four form factors we write 



*" = z + ik F - (46) 



and note that the cross section is a function of six form factors. 
Then the one-loop results apparently equals: 



[ Z = 1] _ * [Z = 0]. (4.7) 



d cos -d d cos -d d cos $ 

5 QED annex 

5.1 QED vertices and soft photon contributions 

Here we present virtual corrections due to QED vertices, a factorised part due to QED boxes 
and soft photon contributions. The expressions of this subsection can be also casted from [8] . 
The formal structure of factorised virtual and soft contributions is as follows: 



irt + soft 



*o2 Cvirt + soft I r\ r\ Cvirt+soft i /^)2 Cvirt+soft 

7T 



(5.1) 



There are three types of contributions: ISR, FSR and IFI. 
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5.1.1 Initial state radiation (ISR) 

Contributions of the initial state QED e + e~7 vertex and ISR soft are short, since electron 
mass is ignored: 



rvirt 
°ISR 



srsoft 



mz 



1 



2 3^ 
1 



-ln^f \l e - 1 ) - -l 2 e + -l e - 2 + 4Li 2 (1) , 



(5.2) 



where 



k = In — 

\ m e, 



(5.3) 



5.1.2 Initial— final state interference (IFI) 

This originates from contributions of QED boxes: 77, Z7 and initial-final state soft photons 
interference: 



with 



^soft 



S,t) 



s virt 

"IFI 



STSoft 



2 e 



-2 In — In — . 

A 2 W_ 

21n — In — 

A 2 -u_ 



F soft (s,t) -F sott (s,u) 



soft 



- In 2 77 + 2 In 77 In 1 + 

2 y p+t- 



-In 2 (l + 
+2Li 2 (1 - 



/M- 
2* 



2Li 2 



-/3_t_ 



(5.4) 
(5.5) 

(5.6) 



/3+t_ + 2m 2 



2Li 2 f-M -2Li 2 (l), 



where we introduce the notations 



Am? 



/?+ = ! + /?, P- = l-P, 



rj 



(5.7) 



5.1.3 Final state radiation (FSR) 

Contributions of one- loop QED ffj vertex and final state soft photon radiation are: 



"FSR 



1 m t 
- ln A^ 

2/3 



:i+i 2 

2/5 



In rj — 1 



- -/3 In 77 -2 



- - In 2 r] + 2 In 77 In (1 - rj) + 2Li 2 (77) + 4Li 2 (1) 
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2 r 



rsoft 



In 

+■ 



4oo 
(1 + /^ 2 ) 



(1 + /3 2 ) 
2/5 



In 77 — 1 



2/9 



- - In 2 77 + 2 In 77 In (1 - 77) + 2Li 2 (77) - 2Li 2 (1) 



(5.8) 



Contribution of the ISR Eq. (5.2) may be received from these expressions in the limit m t = 
m P — > 0. 



5.1.4 Non-factorized final state vertex 'anomalous' contributions 

For presentation of this contribution let us introduce the definition 

(3-1 



In 



13 + 1' 



(5.9) 



The 'anomalous' part of QED vertex contribution to the differential cross-section reads: 



da a 
dcosd 



^N C ^Q 2 



QlQ 2 t + 2Q e Q t v e v t Re( Xz ) + (v 2 e + a 2 ) v 2 \xz\ 2 ) (st + t 2 _) Re{L n ) 



+Q e Q t a e a t s (s + 2£_) Re(L n x z ) 



+ ((,>! + 



2 \ 2 
a A a t 



(s - 4m 2 ) + 2 (st + t 2 .) 



+2?j e a e ?j t a t s (s + 2t_) )\x z \ 2 Re(L n ) 



5.2 An alternative form of the cross-section for QED boxes 



(5.10) 



Here we present some useful formulae which are not in the main stream of our approach 
(described in previous Sections), but that were used for internal cross checks of calculations 
of the QED part of the process under consideration. 

The QED boxes Eqs. (2.26) and (2.27) may be greatly simplified purely algebraically. For 
the sum of AA and ZA boxes one may easily derive the cross-section: 

rfcr BOX 2a 3 



dcosfi 



-(3Q e Q t N c Re \ Q 2 e Q 2 T v + Q e Q t Xl 



v e v t (J 7 * +H V ) + a e a t (JF* + Q v ) 



+ \Xz\ 2 {yl + a 2 ) {v 2 H v + a 2 H A ^ + 2a e v e a t v t (Gv + Ga) }, 
where Xz{s) is defined by Eq. (1.4) and the six cross-section form fartors are: 



(5.11) 







(«), 


T A 


= T A {t)+T A 


(«), 


T~C V 


= TL V (t) — 7i v 


(«), 


H A 


= H A (t) — H A 


(«), 


Gv 


= Gv(t) + Gv( 


u), 


G A 


= G A (t) + G A ( 


u) , 



(5.12) 
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with 



n A (t) 



-{- 

si 4 



2m 2 + (s + 2t_) 

2 



smf 



2 

1 f s + 2i_ 



(s + t_) 



^aa(-s, -t;m e ,m t ) 

sC ( — ml, —ml, —s; 0, m e , 0) + (s — 4m 2 ) C ( — m 2 , — m 2 , — s; 0, m t , 0) 

2t72 2 r ] \ 

— ± 2m 2 t Co(-m 2 t ,-m 2 t ,-s;0,m t ,0) + <(-m 2 ;m t ,0) - <( - s; 0, 0) 

ZA 3r L J / 

Bjf ( - t; m e , m*) - 5 F ( - m 2 ; m t , 0) 
E F (-*;m e ,m t ) -Bft-s; 0,0)] J, 

~ *— -*; m e , "it) - m 2 QsC ( - m\ , -m\, -s; 0, m t , 0) 
) 5 F ( - t; m e , m t ) - -Bjf( - m 2 t ; m t , 0) ^ 
<(-t;m e ,m t )-<(- S ;0,0)]}, 

lr r l 

W o (0 = -^j - {f- + (s + t-) 2 )t- Jza(s, -t;m e ,m t ) - -J AA (-s, -t;m e ,m t ) 
+C { - mf, -ml, -t; m t , M z , m e ) 

+s- (t- (s+ + 2t) ^J ZA (-s, -t; m e , m t ) + C ( - m 2 , -m 2 , -t; m t , M z , m e ) 



(5.13) 



(5.14) 



M 2 C ( - m 2 , -m 2 , -t; m t , M z , m e ) - 2Bq( - t; m e , m t ) 
+B%( - m 2 t ; m t , 0) + B^{ - m\; M z , m t ) 



-st 



H v (t) = H (t) + 



C ( - m 2 , -m 2 , -s; M z , m e , 0) + C ( - m 2 , -m 2 , -s; M z , m t , 0) 
<(-t;m e ,m t )-<(-s;M z ,0) )}, 

Jza{s, -t; m e , m t ) - ^ J AA (-s, -t; m e , m t ) 



s + *_) 

2m? 



(5.15) 



+C ( - mf, -m 2 e , -t; m t , M z , m e ) 



+s-t( - C ( - m\,-m\,-s; M z ,m t ,0) + —— s_C ( - ml, -ml, -s; M z ,m t ,0) 
\ A 3r L 

+B F ( - m 2 t ; m t , 0) + 5 F ( - m 2 ; M z , m t ) - 2B F ( - s; M z , 0)] ) }, (5.16) 
M 2 J ZA (-s, -t; m e ,m t ) 



„, , s 2m? I 
Wo(*) + -pM* 
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s (^Jaa(s, -t; m e , m t ) - C ( - m 2 , -m 2 e , -t; m t , M z , m e )^j 

+s_ - tC ( - m 2 , -m 2 , -s; M z , m t , 0) + s + C ( - m 2 , -m 2 e , -t; m t , M z , 

-B^( - t; m e , m t ) + Bg( - s; M z , 0) J, 



(5.17) 



G v (t) = — j(s + 2t_)t_ 



^ J A a{s, -t; m e , m t ) + C ( - m 2 , -m 2 e , -t; m t , M z , 



Jza(s, -t;m e ,m t ) 



(s+ + 2t_) ^J ZA {-s, -t; m e , m t ) + C ( - m 2 , -m 2 e , -t; m t , M z , m e ) 

C ( - ml, -ml, -s; M z , m e , 0) + C ( - m 2 t , -m 2 t , -s; 

1 r „ 



M z ,m t ,0) 



+ \ M K 

+m 2 t ^2C ( - m\, -m 2 e , -t; m t , M z , m e ) + — M 2 Z C Q ( - m\, -m 2 e , -t; m t , M z , 
3q( - t; m e , m t ) + B${ - m\; m t , 0) + B%( - m 2 ; M z , m t ) 
' 7 ' ) [ <( - /: m e , m t ) -B^(-s;M z , 0) 



m e ) 



-2B* 



(5.18) 



G A (t) = Q v {t)-m 2 t S — 

s z 



s z L 

+s+C ( - ra t 2 , -ra 2 , -s; M z , m t , 
+2B^(-t;m e ,m t )-B^ 



, z ,-m z e ,-t;m t ,M z ,m e ) 



0) 



? F ( - m 2 ; M z , rm) - <( - m 2 ; m t , 0) 



(5.19) 



The Eqs. (5. 11)— (5. 19) were coded as a separate branch of eettLib and together with vertex 
QED contributions described in the previous subsections was used for internal cross-check of 
QED part of calculations. 

Some factorized part of the AA and ZA boxes contribution is not included in Eq. (5.11). 
It has a form 



da a 
dcost? 71 



(5.20) 



where 5™* is given by Eq. (5.4). 

The whole QED contribution can be written as follows 

^ QED _ ^ BOR Y irt+ soft da a da™ 



+ 



+ 



dcosfi dcosfi dcosi!) dcosfi' 

where <5 virt + soft is defined by Eqs. (5.1)-(5.8). 



(5.21) 
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6 Numerical results and discussion 



All the formulae derived in this paper as well as in Ref. [3] are realized in a FORTRAN code with 

a tentative name eef f Lib. Numbers presented in this section are produced with updated, 

February 2002, version of the code. As compared to December 2000 version, used to produce 

numbers for Ref. [3], current version contains full QED corrections together with the soft 

photon contribution to the angular distribution da / dcos'd. Morever, two bugs of December 

2000 version were fixed, which resulted in a change of numbers. First, for light final state 

fermion masses, the numerical precision was being lost. After curing this oddity, the agreement 

between eef f Lib and ZFITTER numbers became even better. Secondly, there was a bug in 

a part of FORTRAN code computing ZZ box contribution. Its fixing resulted in a change of 

numbers (in 4th-5th digits) for the case of heavy final state fermion masses (top quark). 

Since numbers which were presented in Ref. [3] are changed anyway, we decided to present 

again all the Tables that were already given in Ref. [3]. On top of it, we will show several 

new examples of numbers. In particular, we will show a comparison of the electroweak form 

factors (EWFF) including QED corrections between eeffLib and another FORTRAN code, 

which was automatically generated from form log files with the aid of a system s2n_f {symbolic 

to numbers), producing a FORTRAN source code - - a part of our CalcPHEP system. This 

comparison provides a powerful internal cross-check of our numerics that practically excludes 

appearance of bugs of a kind discussed above. 

We begin with showing several examples of comparison with ZFITTER v6.30 [9]. In the 

present realization, eeffLib does not calculate M w from \i decay and does not precompute 

either Sirlin's parameter Ar or total Z width, which enters the Z boson propagator. For 

this reason, the three parameters: M w , Ar , V z were being taken from ZFITTER and used as 

INPUT for eeffLib. Moreover, present eeffLib is a purely one- loop code, while in ZFITTER 

it was not foreseen to access just one- loop form factors with users flags. To accomplish the 

goals of comparison at the one-loop level, we had to modify the DIZET electroweak library. 

The most important change was an addition to the SUBROUTINE R0KANC: 
* 

* For eett 
* 

FLL=(XR0K(1)-1D0+DR )*R1/AL4PI 
FQL=FLL+ (XR0K (2) -1D0) *R1/AL4PI 
FLQ=FLL+ (XR0K (3) -1D0) *R1/AL4PI 
FQQ=FLL+ (XR0K (4) -1D0) *R1/AL4PI 

with the aid of which we reconstruct four form factors from ZFITTER's effective couplings p 
and k's (F ld and F QD do not contribute in massless approximation). 

6.1 Flags of eeffLib 

Here we give a description of flags (user options) of eeffLib. While creating the code, we 
followed the principle to preserve as much as possible the meaning of flags as described in the 
ZFITTER description [10]. In the list below, a comment 'as in ZFD' means that the flag has 
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exactly the same meaning as in [10]. Here we describe an extended set of flags of February 
2002 version of eef fLib. 



ALEM=3 
ALE2=3 
VP0L=0 



as in ZFD 
as in ZFD 

=0 a(0); = 



1,=2 as in ZFD; =3 is reserved for later use 



Note that the flag is extended to VP0L=0 to allow calculations 'without running of a'. 



• 


QCDC=0 


as 


in 


ZFD 


• 


IT0P=1 


as 


in 


DIZET (internal flag) 


• 


GAMS=1 


as 


in 


ZFD 


• 


WEAK=1 


as 


in 


ZFD (use WEAK=2 in v6.30 to throw away some higher order terms) 


• 


IM0MS=1 


=0 


a- 


scheme; =1 GFermi-scheme 



New meaning of an old flag: switches between two renormalization schemes; 
B0XD=6 

Together with WEAK=0 is used for an internal comparison of separate boxes and QED 
contributions: 

B0XD ! =1 with 77 boxes 
=2 with Z7 boxes 
=3 with 77 and Z7 boxes 
=4 with all QED contributions 
Together with WEAK=1 (working option), it has somewhat different meaning: 
B0XD ! =0 without any boxes 
=1 with 77 boxes 
=2 with Z7 boxes 
=3 with 77 and Z7 boxes 
=4 with WW boxes 
=5 with WW and ZZ boxes 
=6 with all QED and EW boxes 

'Treatment' options. 

GAMZTR=1 treatment of T z . 

The option is implemented for the sake of comparison with FeynArts: 
GAMZTR=0 T z = 
GAMZTR=1 T z ^ 

EWFFTR=0 treatment of EW form factors. 

Switches between form factors and effective ZFITTER couplings p and re's. The option is 
implemented for comparison with ZFITTER: 
EWFFTR=0 electroweak form factors 
EWFFTR=1 effective coullings p and re 
FERMTR=1 treatment of fermionic masses. 

Switches between three different sets of 'effective quark masses': 
FERMTR=1 a 'standard' set of fermions masses 
FERMTR=2 , 3 'modified' sets 
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• VP0LTR=1 treatment of photonic vacuum polarization. 

Switches between lowest order expression, a(s) = a [1 + Aa(s)], and its 'resummed' 
version, a(s) — a/ [1 — Aa(s)]: 
VP0LTR=0 lowest order 
VP0LTR=1 resummed 

• EWRCTR=2 treatment of electroweak radiative corrections. 
Switches between three variants for vertex corrections: 
EWRCTR=0 electroweak form factors contain only QED additions 
EWRCTR=1 electroweak form factors do not contain QED additions 
EWRCTR=2 electroweak form factors contain both QED and EW additions 

• EMASTR=0 treatment of terms with Yti(s/m 2 e ) in 77 and Z7 boxes, which are present in 
various functions but cancel in sum: 

EMASTR=0 these terms are suppressed in all functions which they enter 

EMASTR=1 these terms are retained in all functions which results in loosing of computer 

precision owing to numerical cancellation; results for EMASTR=0 and EMASTR=1 are equal 

• EWWFFV=1 treatment of vertex and box diagrams with virtual W boson, switches between 
two variants: 

EWWFFV=0 variant of formulae without 6-quark mass 
EWWFFV=1 variant of formulae with finite 6-quark mass 

Options affecting QED contributions. 

• IQED=4 variants of inclusion virtual and soft photon QED contributions: 
IQED=1 only initial state radiation (ISR) 

IQED=2 only initial-final interference (IFI) 

IQED=3 only final state radiation (FSR) 

IQED=4 all QED contributions are included 

• IB0X=4 is active only if IQED=2 or 4 and affects only Eq. (5.11): 
IB0X=0 AA boxes interfering with 7 exchange BORN 

IB0X=1 AA boxes 
IB0X=2 ZA boxes 
IB0X=3 or 4 AA+ZA boxes 

6.2 eef f Lib— ZFITTER comparison of scalar form factors 

First of all we discuss the results of a computation of the four scalar form factors, 

F LL (s,t), F QL (s,t), F LQ (s,t), F QQ (s,t), (6.1) 

for three variants: 

1) without EW boxes, i.e. without gauge-invariant contribution of ZZ boxes, and without 
f = 1 part of the WW box, Eq. (1.2.93); 

2) with WW boxes; 

3) with WW and ZZ boxes. 
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Table 1: EWFF for the process e+e~ -> uu. eef f Lib-ZFITTER comparison. 



Without EW boxes 




100 GeV 


200 GeV 


300 GeV 


FF 


i ^ 










M w /10 
M w 
10M W 


13.47773 - ^1.84784 
loAl 1(6 — 2l.o4/o4 
13.47773 - il.84784 


16.22034 -?10.49412 

1 (? oonQ/i ;m /i O/f i o 
10.ZZUo4 — ZlU.4y41z 

16.22034 -il0.49412 


23.75240 -ill. 27469 

no vcro/in ;i 1 07/1(?n 

23.75240 -ill.27469 


ZFITTER 


13.47771 - il.84786 


16.22031 -ilO.49405 


23.75237 -ill.27464 


Fql 


M w /10 
M w 
10M W 


29.34721 + ^3.67330 

on q^oi i ao C700c\ 
Z\).o4:(Zl + 16. 6/330 

29.34721 + i3.67330 


30.33891 + ^3.34531 
30.33oyl + 16. 34531 
30.33891 + i3.34531 


31.64553 + «2. 75258 

oi a tr tr o i ^'o ^7txotxo 
ol. 64553 + 2Z./5z5o 

31.64553 + il. 75258 


ZFITTER 


29.34720 + i3.67330 


30.33889 + i3.34535 


31.64552 + i2. 75259 




M w /10 
M w 
10M W 


29.13302 + ?3. 26972 

on i oono i ao nen^o 
zy.looOz + 16.ZW/Z 

29.13302 + i3.26972 


30.03854 + il. 54158 

30.03o54 + «1.04l0o 

30.03854 + il.54158 


31.68636 - lO. 22635 

01 dQCiOCi ACl OO^OtX 

31.6o636 — i\).ZZb6o 
31.68636 - i0.22635 


ZFITTER 


29.13304 + i3.26973 


30.03855 + il.54163 


31.68635 - i0.22634 


Fqq 


M w /10 
M w 
10M W 


44.90390 + i8.85688 
44.90390 + i8.85688 
44.90390 + i8.85688 


43.80286 +il0.02412 
43.80286 +il0.02412 
43.80286 +^10.02412 


44.21223 +il0.83899 
44.21223 +il0.83899 
44.21223 +^10.83899 


ZFITTER 


44.90392 + i8.85688 


43.80285 + il0.02411 


44.21224 +il0.83894 


WW is added 


F LL 


M w /10 
M w 
10M W 


12.94469 - il.84784 
12.94469 - il.84784 
12.94469 - il.84784 


9.34066 - i9.42482 
9.34066 - i9.42482 
9.34066 - i9.42482 


9.03908 — ill.55971 
9.03908 — zll.55971 
9.03908 — zll.55971 


ZFITTER 


12.94468 - il.84786 


9.34065 - - i9.42467 


9.03903 — zll.55958 



In this comparison we use flags as in subsection 6.1 and, moreover, 

M w = 80.4514958 GeV, 
Ar = 0.0284190602, 

Y z = 2.499776 GeV. (6.2) 

In Table 1 we show an example of comparison of four form factors F LL QL LQ QQ (s, t) between 
the eeffLib, where we set m t = 0.2 GeV and ZFITTER (the latter is able to deliver only 
massless results). The form factors are shown as complex numbers for the three c.m.s. energies 
(for t = — s/2) and for the three values of scale \i = M w /10, M w , 10M W . The table 
demonstrates scale independence and very good agreement with ZFITTER results (6 or 7 digits). 
One should stress that total agreement with ZFITTER is not expected because in the eeffLib 
code we use massive expressions to compute the nearly massless case. Certain numerical 
cancellations leading to losing some numerical precision are expected. We should conclude 
that the agreement is very good and uniquely demonstrates that our formulae have the correct 
m t — * limit. 
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Table 2: EWFF for the process e + e~ -> mm. eef f Lib-ZFITTER comparison. 



With ZZ boxes 




100 GeV 


200 GeV 


300 GeV 


FF 


u 










M w /10 
M w 
10M W 


12.89584 - il.84784 
12.89584 - il.84784 
12.89584 - il.84784 


8.24736 -ilO.64666 
8.24736 -ilO.64666 
8.24736 -ilO.64666 


8.98375 -il2.88478 
8.98375 -il2.88478 
8.98375 -il2.88478 


ZFITTER 


12.89583 - il.84786 


8.24736 -il0.64651 


8.98370 -il2.88466 


Fql 


M w /10 
M w 
10M W 


29.30447 + i3.67330 
29.30447 + i3.67330 
29.30447 + i3.67330 


29.38219 + i2.27610 
29.38219 + i2.27610 
29.38219 + i2.27610 


31.59711 + il.59302 
31.59711 + il.59302 
31.59711 + il.59302 


ZFITTER 


29.30445 + i3.67330 


29.38216 + i2.27613 


31.59710 + il.59304 


Flq 


M w /10 
M w 
10M W 


29.10829 + i3.26972 
29.10829 + i3.26972 
29.10829 + i3.26972 


29.48510 + i0.92307 
29.48510 + i0.92307 
29.48510 + i0.92307 


31.65836 - - i0.89713 
31.65836 - i0.89713 
31.65836 - i0.89713 


ZFITTER 


29.10832 + i3.26973 


29.48512 + i0.92312 


31.65835 -- i0.89711 


Fqq 


M w /10 
M w 
10M W 


44.88226 + i8.85688 
44.88226 + i8.85688 
44.88226 + i8.85688 


43.31855 + i9.48287 
43.31855 + i9.48287 
43.31855 + i9.48287 


44.18773 +il0.25200 
44.18773 +il0.25200 
44.18773 +il0.25200 


ZFITTER 


44.88228 + i8.85688 


43.31854 + i9.48286 


44.18773 +il0.25196 



In Table 2 we show a similar comparison with ZFITTER when ZZ boxes are added. As 
seen, the agreement has not deteriorated. 

6.3 eef f Lib— ZFITTER comparison of IB A cross-section 

As the next step of the comparison of eef f Lib with calculations from the literature, we present 
a comparison of the IBA cross-section. 

In Table 3 we show the differential cross-section Eq. (4.44) in pb for three values of cos-# = 
—0.9, 0, +0.9, with input parameters of Eq. (6.2) and with constant e.m. coupling a = a(0). 

Next, we present the same comparison as in Table 3, but now with running e.m. coupling. 
Since the flags setting VP0L=1, which is relevant to this case, affects ZFITTER numbers, we 
now use, instead of Eq. (6.2), the new INPUT set: 

M w = 80.4467671 GeV, 
Ar = 0.0284495385, 

Y z = 2.499538 GeV. (6.3) 

The numbers, shown in first two rows of Tabs. 3 and 4 exhibit a very good level of agreement 
for light quark masses, while the third rows illustrate the mass effect due to heavy top. 

Finally, in Table 5, we give a comparison of the cross-section integrated within the angular 
interval | cos^| < 0.999. (Flags setting is the same as for Table 4.) 
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Table 3: IBA, First row - ZFITTER (uu channel); second row - eeffLib (mt = 0.1 GeV, 
mb = GeV); third row - eeffLib (m t = 173.8 GeV); with constant e.m. coupling a = a(0). 





100 GeV 


200 GeV 


300 GeV 


400 GeV 


700 GeV 


1000 GeV 


costf = -0.9 


4 1 .00400Z 


n oni qoq 
U.zylozo 


u.ioyoiu 








47.661843 


0.291827 


0.169515 


0.103284 


0.035319 


0.017204 








0.162193 


0.044088 


0.018927 


cos = 


59.768387 


1.718830 


0.695061 








59.770299 


1.718870 


0.695075 


0.376871 


0.117279 


0.055873 








0.264713 


0.112918 


0.054209 


costf = 0.9 


168.981978 


5.954048 


2.292260 








168.991144 


5.954182 


2.292302 


1.222354 


0.372912 


0.176038 








0.438453 


0.293399 


0.154785 



Table 4: IBA, First row - ZFITTER (uu channel); second row - eeffLib (mt = 0.1 GeV, 
nib = GeV); third row - eeffLib (m t = 173.8 GeV); with running e.m. coupling a = a(s). 





100 GeV 


200 GeV 


300 GeV 


400 GeV 


700 GeV 


1000 GeV 


cos $ = -0.9 


45.404742 
45.404598 


0.386966 
0.386966 


0.225923 
0.225923 


0.138065 
0.194752 


0.048621 
0.058013 


0.024156 
0.025959 


cos $ = 


60.382423 
60.382562 


1.882835 
1.882837 


0.771939 
0.771939 


0.421410 
0.303683 


0.133475 
0.130173 


0.064245 
0.062838 


costf = 0.9 


173.467517 
173.467543 


6.450000 
6.450000 


2.510881 
2.510881 


1.346620 
0.492546 


0.417295 
0.330401 


0.198842 
0.175564 



A typical deviation between eeffLib and ZFITTER is of the order ~ 10~ 6 , i.e. of the 
order of the required precision of the numerical integration over cos-$. Examples of numbers 
obtained with eeffLib, which were shown in this section, demonstrate that ZFITTER numbers 
are recovered for light m t . 

We conclude this subsection with a comment about technical precision of our calculations. 
We do not use looptools package [4]. For all PV functions, but one, namely D function, we 
use our own coding where we can control precision internally and, typically, we can guarantee 
11 digits precision. For Dq function we use, instead, REAL* 16 T0PAZ0 coding [11] and the only 
accessible for us way to control the precision is to compare results with those computed with 
looptools package. This was done for a typical D functions entering ZZ box contributions. 
We got an agreement within 14-15 digits between these two versions for all = 400 — 10000 
GeV and cosv* = 0.99, 0, -0.99. 
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Table 5: eef f Lib(L)-ZFITTER(Z) comparison of the total cross-section. Cross-sections are 
given in pb: the first row - a^ ot , i.e. eef f Lib (m t = 0.1 GeV); the second row - a^ t , i.e. 
ZFITTER (uu channel); the third row shows the absolute deviation a^ ot — af ot . 



100 GeV 


200 GeV 


300 GeV 


Ctot 


0"fb 


Ctot 


0"fb 


Ctot 


0"fb 


160.8981 


70.98416 


5.021810 


3.360848 


2.031754 


1.269556 


160.8980 


70.98406 


5.021808 


3.360848 


2.031754 


1.269556 


0.0001 


0.00010 


0.000002 


0.0 


0.0 _j 


0.0 ± 



6.4 Comparison with a code generated by s2n_f 

Here we present a numerical comparison of the complete scalar form factors Eq. (3.1) extracted 
from two independently created codes: 'manually written' eettLib and a code, 'automatically 
generated' by s2n_f software. We use a special input parameter set here: all lepton masses a 
and a conversion factor from GeV -1 to pb are taken from 2000 of Particle Data Tables while 
for quark and photon and gauge boson masses we use: 

m uAc ^ h = 0.062, 0.083, 1.50, 0.215, 173.8, 4.70 GeV, 

A = 1 GeV, M z = zm = 91.1867 GeV, M w = 80.4514958 GeV. (6.4) 



Table 6: EWFF for the process e + e — > tt. eettLib, first rows; s2n_f , second rows. 



v/5 



400 GeV 



700 GoV 



1000 GoV 



Fll 



Flq 
F QQ 
Fld 
Fqd 



68.36399900074 
68.36399900068 
75.12465846647 
75.12465846641 
81.01546270426 
81.01546270420 
225.63977621858 
225.63977621832 
-0.57522852857 
-0.57522852857 
0.16677424366 
0.16677424366 



1.24743850729 
1.24743850728 
+ i34.81991916400 
+ i34.81991916400 
+ il9.81343626967 
+ il9.81343626968 
+ il54.37838168488 
+ il54.37838168491 
+ iO.34010611241 
+ iO. 34010611241 
- iO. 34326069364 
- iO. 34326069364 



79.63957322115 
79.63957322113 
76.19283172015 
76.19283172013 
82.67283873006 
82.67283873004 
207.09189805263 
207.09189805254 
-0.33030593699 
-0.33030593699 
0.29925308488 
0.29925308488 



- i20. 53758995637 

- i20.53758995637 
+ 128.44336684106 
+ i28.44336684106 
+ 113.79952080171 
+ 113.79952080171 
+ 1133.45188150116 
+ 1133.45188150117 
+ iO. 14897150833 

+ iO. 14897150833 
- iO. 14107543098 

- iO. 14107543098 



80.47816819240 - 
80.47816819239 ■ 
75.95332822621 - 
75.95332822620 - 
83.26485989744 - 
83.26485989743 - 
194.07155316803 - 
194.07155316799 - 
-0.22418674728 ■ 
-0.22418674728 ■ 
0.23436559470 ■ 
0.23436559470 ■ 



i26.71016937527 
126.71016937527 
i27.77201429453 
i27.77201429453 
il2.23741074712 
il2.23741074712 
1134.33226297675 
1134.33226297675 
iO. 08847119487 
iO. 08847119487 
iO. 05839137636 
iO. 05839137636 



Fll 
Fql 
Flq 
Fqq 
Fld 
Fqd 



48.42950001713 
48.42950001707 
68.02678564355 
68.02678564349 
73.37133716227 
73.37133716220 
196.60425612149 
196.60425612123 
-0.56319765502 
-0.56319765502 
0.15893936555 
0.15893936555 



+ 18.26103890366 
+ i8.26103890367 
+ 137.08805801477 
+ i37.08805801477 
+ i22.69397728402 
+ i22.69397728403 
+ il62.74818773960 
+ il62.74818773963 
+ iO. 33645326768 
+ iO. 33645326768 

- iO. 37254018572 

- iO. 37254018572 



28.23570422021 
28.23570422019 
58.00469565609 
58.00469565607 
62.40775508619 
62.40775508616 
132.63279537966 
132.63279537957 
-0.29067043403 
-0.29067043403 
0.26429138671 
0.26429138671 



+ i2.43705570966 
+ i2.43705570966 
+ i33.82433896562 
+ i33.82433896562 
+ 120.75544388763 
+ 120.75544388764 
+ 1152.68259938740 
+ 1152.68259938741 
+ iO. 13992893252 
+ iO. 13992893252 
- iO. 15437851127 
- iO. 15437851127 



16.54896558498 - 
16.54896558497 ■ 
52.56343218854 H 
52.56343218853 H 
56.94788960099 H 
56.94788960098 H 
98.94491326876 4 
98.94491326872 -\ 
-0.18096486789 - 
-0.18096486789 - 
0.18981891199 - 
0.18981891199 - 



- iO. 77082434583 

- iO. 77082434583 
i34.28418004972 
i34.28418004972 
i20. 65389886886 
i20. 65389886886 
il57.45863002555 
il57.45863002556 
iO. 08187546112 
iO.08187546112 
iO. 06181088399 
iO. 06181088399 



Fll 
Fql 
Flq 
Fqq 
Fld 
Fqd 



35.17736865724 
35.17736865718 
61.03099608330 
61.03099608324 
66.08215572935 
66.08215572929 
167.63393504156 
167.63393504130 
-0.56772633347 
-0.56772633347 
0.18031346246 
0.18031346246 



+ il4.84038724783 
+ 114.84038724784 
+ i39.09196533610 
+ i39.09196533611 
+ i25.04151178684 
+ i25.04151178685 
+ il70.36384103672 
+ il70.36384103675 
+ iO. 34299744419 
+ iO. 34299744419 

- iO. 40091423652 

- iO. 40091423652 



0.21531292996 
0.21531292994 
40.77942026097 
40.77942026095 
44.50915974057 
44.50915974055 
59.87568281297 
59.87568281288 
-0.32035310873 
-0.32035310873 
0.34968026058 
0.34968026058 



+ il3.66645015866 
+ 113.66645015866 
+ i37.94118444135 
+ i37.94118444136 
+ 125.51875704261 
+ i25.51875704261 
+ il68. 13599380718 
+ il68. 13599380719 
+ iO. 14419510235 
+ iO. 14419510235 
- iO. 16945266925 
- iO. 16945266925 



-18.22896674792 ■ 
-18.22896674793 ■ 
30.73687048410 - 
30.73687048409 - 
34.09235431695 - 
34.09235431694 - 
7.10290370391 A 
7.10290370387 4 
-0.21547870582 - 
-0.21547870582 - 
0.29109057806 - 
0.29109057806 - 



- il:S. 10520552155 

- 113.10520552155 
i39. 06896169999 
i39. 06896169999 
i26.17012224125 
i26.17012224125 
il75.24101109592 
il75. 24101109593 

■ iO. 08254457292 

■ iO. 08254457292 
iO. 06775284666 
iO. 06775284666 
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As seen from the table numbers agree within 11-13 digits, i.e. REAL*8 computational 
precision is saturated. Form factors F LDiQD are multiplied by 10 4 to make more digits visible. 

Next Table demonstrates eettLib-s2n_f comparison for the complete one-loop differential 
cross-sections da^ 1 )/d cos § for the standard input set. As seen, numbers agree within 12-13 
digits. 



Table 7: 4^ — * for the process e + e — > tt. eettLib-s2n_f 

a cos 17 ^ 





400.0 


700.0 


1000.0 


cosi/ 1 








-0.900 


0.22357662754774 
0.22357662754769 


0.06610825350063 
0.06610825350063 


0.02926006442715 
0.02926006442715 


0.000 


0.34494634728716 
0.34494634728707 


0.14342802645636 
0.14342802645634 


0.06752160108814 
0.06752160108813 


0.900 


0.54806778978208 
0.54806778978194 


0.33837133344667 
0.33837133344664 


0.16973989931024 
0.16973989931023 



6.5 About a comparison with the other codes 

As is well known, the one-loop differential cross-section of e + e~ — > tt may be generated with 
the aid of the FeynArts system [4]. Previous attempt to compare with FeynArts are described 
in [3]. In December 2001, we were provided with the numbers computed with the FeynArts 
system [13] for da/dcostf with and without QED contributions at y/s = 700 GeV and three 
values of cosfl = 0.9, ,0, —0.9. After debugging of our code eettLib, as described in the 
beginning of this section, we eventually reached 11 digits agreement both for the tree level 
and one-loop corrected cross-sections. 

We do not update Fig. 13 and Fig. 14 of [3], since the differences with updated version is 
not seen. 

Table 8: \ for the process e+e~ — > tt with soft photons, £™ ax = y/s/10. 





400.0 


700.0 


1000.0 


COS1/ 1 








-0.900 


0.17613018248935 


0.05199100267864 


0.02310170508071 


-0.500 


0.21014509428358 


0.06560630503586 


0.02882301902010 


0.000 


0.27268108572063 


0.11496514450150 


0.05495088904853 


0.500 


0.35592722356682 


0.19615154401629 


0.09941700898317 


0.900 


0.43637377538440 


0.27915043976042 


0.14426233253975 



Recently, a Bielefeld-Zeuthen team [5] performed an alternative calculations using the 
DIANA system [12]. Working in close contact with this team, we managed to perform several 
high-precision comparisons reaching for separate contributions an agreement in 10 digits. 
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The results of a comparison between FeynArts and Bielefeld-Zeuthen team are presented 
in detail in [6]. 

As another example we present in Table 8 the same cross-section [da / d cos $] SM as given in 
tables of [6]. For the complete cross-section, including soft photons, we agree with Bielefeld- 
Zeuthen calculations within 7-8 digits. 
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